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Abstract We propose to use nonparametric Bernstein copulas as bivariate pair-copulas in high- 
dimensional vine models. The resulting smooth and nonparametric vine copulas completely obvi- 
ate the error-prone need for choosing the pair-copulas from parametric copula families. By means 
of a simulation study and an empirical analysis of financial market data, we show that our pro- 
posed smooth nonparametric vine copula model is superior to competing parametric vine models 
calibrated via Akaike's Information Criterion. 
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1 Introduction 



Following the growing criticism of elliptical models, copulas have emerged both in insurance 
and risk management as a powerful alternative for mod eling the comple t e dep endence structure 



Embrechts et al. 



(l2002h . the literature on 



of a multivariate distribution. Since the seminal work by 
copulas and their use in risk management applications has grown exponentially w ith several s tudies 
concentrating on sta tistical inference and model sele ction for copulas (see, e.g., 



Genest et al. 



2012 



Yea et alibom fl 



2009bh a s well as applications (see, e.g. JChan and Kroesel . l2012l : lGrundke and PoUe . 



Kim et al, 



2007 



As simple parametric mq d .els are often not flexible enough t o model hig h-dir nensional distri 



Whelan 



( 2004 



butions, recent works by Hod (I1996L 119971) . Bedford and Cookd (1200 IL 120021) and 
have proposed copula models which are highly flexible but at the same time still tractable even 
in higher dimensions. Most notably, vine copulas (also called pair-copula constructions, PCC in 
short) have emerged as the most promising tool for modeling dependence structures in high dimen- 
sions]^ Vine copulas consist of a cascade of conditional bivariate copulas (so called pair-copulas) 
which can each be chosen from a different parametric copula family. As a result, vine copulas 



are extremely flexible yet still tractable even in high dimens ions as al 



Aasetal 



computations necessary in 



2009, for a first discussion 



statistical inference are performed on bivariate data sets (see 
of vine copulas in an applied setting). 

Similar to the bivariate caseB the correct selection of the parametric constituents of the vine, 
i.e., the pair-copulas, is crucial for the correct specification of a vine copula model. In case of 
the popular C- and D-vine specifications, the calibration and estimation of a (i-dimensional vine 
requires the selection and estimation of d{d — l)/2 different pair-copulas from the set of candidate 
bivariate parametric copula families. Thus, a vine model's increased flexibility only comes at the 
expense of an increased model risk. 

As a remedy, recent studies have suggested to select the param etric pair-copulas based on 
graphical data inspection and goodness-of-fit tests (lAas et al.l |2009|) a nd to employ sequential 



heuristics based on Akaike' s Information Criter ion ( AIC) (see, e.g.. iBrechmann et al.L 12012 : 



DiBmann et al, 



forthcoming!) . iKurowickal (|2010r) and 



Brechmann et al 



(|2012h propose strate- 



gies for simplifying vines by replacing certain pair-copulas by the independence copula (yield- 
ing a truncated vine copula) or the Gaussian copula (yielding a simplified vine). Finally, 



'a literature review with a special emphasis on finance-related papers using copulas is given by iGenest et al 
(l2009ah . An overview of the different branches of the copula literature is given by Embrechts (2009) 



^C ompeting modeling concepts l ike nested and hierarchical Archimedean copulas are analyzed by Aas and Berd 
( 120091) as well as lFischer et al. I (l2009h . They conjecture that vine copulas should be preferred over nested or hierarchi- 
cal Archi medean copu l as. 

^See lGenest et all (l2009bh . iKole et alJ (l2007h and lWeifil (l201 l[l2012h for discussions of the problem of selecting 
the best fitting parametric copula. 
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Hobaek-Haff and SegersI (|2012h propose the use of empirical pair-copulas in vine models to cir- 
cumvent the problem of selecting parametric pair-copulas. 

In this paper, we u s e the rece n tly p r oposed nonparam etric Bernstein copulas 



(jSancetta and Satchell . 



2004 



Pfeifer et al, 



20091 : biers et all 1201 2h as pair-copulas yielding 



smooth nonparametric vine copula models that do not require the spe c ificati on of parametric 
families. Thus, we extend the ideas laid out by iHobaek-Haff and SegersI (|2012|) by using an ap- 



proximation to the empirical pair-copulas. In contrast to their work, however, we approximate the 
pair-copulas not only nonparametrically but also by the use of continuous functions j In addition . 



espec ially the Bernstein copula has recently attracted attention in insurance modeling (IDiers et al. . 
2012h and has already proven its merits in an applied setting. Therefore, the contributions of the 
proposed smooth nonparametric vine copulas are twofold: First, the use of Bernstein copulas 
completely obviates the need for the error-prone selection of pair-copulas from pre-specified 
sets of parametric copulas. The resulting smooth and nonparametric vine copulas do not only 
constitute extremely flexible tools for modeling high-dimensional dependence structures, they are 
also characterized by a smaller model risk than their parametric counterparts. Second, Bernstein 
copulas have been shown to improve on the estimation of the underlying dependence structure by 
competing nonparametric empirical copulas^The modeling of a vine model's pair-copulas by the 
use of smooth approximating functions is thus a natural extension of recently proposed (highly 
discontinuous) empirical pair-copulas. 

The usefulness of the proposed Bernstein vines is demonstrated by means of a simulation 
study as well as by forecasting and backtesting the Value-at-Risk (VaR) for multivariate portfolios 
of financial assets. 

The results presented in this study show that our proposed vine copula model with smooth non- 
parametric Bernstein pair-copulas outperforms the benchmark model with parametric pair-copulas 
in higher dimensions with respect to the accuracy and numerical stability of the approximation 
to the true underlying dependence structure. While our nonparametric vine copula model yields 
worse average squared errors than a benchmark vine copula calibrated by selecting parametric 
pair-copulas based on AIC values in lower dimensions (e.g., d = 3,5, 7) in our simulations, this 
result is reversed in higher dimensions. For random vectors of dimension d = 11 and higher, the 
parametric benchmark broke down in more than 50% of the simulations due to either the numerical 



^Using smooth functions to approximate the true underlying dependence structure is in line with our intuition. 
The superiori ty of smoo t h non parametric approximations of the copula over simple empirical copulas, however, is 



also found by lShen et all (l2008h . They argue that the improved approximation by the linear B-spline copulas is due to 
their Lipschitz continuity. 

^For example, Bernstein copulas provide a higher rate of consistency than other common nonpara metric estimators 
and do not suffer from boundary bias ( Kulpa , 1999t Sancetta and Satchell , 2004 : Piers et al. , 2012 ). Similarly, other 
approximations as, e.g., linear B-spline copula s have also been s hown to yield lower average squared approximation 
errors than competing discrete approximations (iShen et al.ll2008l) . 
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instability of the parameter and AIC estimation or simply due to extremely inaccurate approxima- 
tions caused by badly selected parametric pair-copulas. In higher dimensions (i.e., the main field 
of application of vine copulas), our nonparametric modeling approach is thus clearly superior to 
a parametric vine copula model. Our risk management application, however, shows that even in 
lower dimensions (d = 5) our nonparametric model yields VaR forecasts that are not rejected 
by a range of formal statistical backtests. Consequently, the slightly worse approximation errors 
of our nonparametric model in lower dimensions do not seem to affect the modeling of a given 
dependence structure too severely thus underlining the usefulness of our proposed model. 

The remainder of this article is structured as follows. Section 2 introduces vine copulas as well 
as Bernstein copulas which we employ as pair-copulas in the vines. Section 3 presents the results 
of a simulation study on the approximation errors of both our nonparametric Bernstein vine copula 
model as well as a heuristically calibrated parametric benchmark model. In Section 4, we conduct 
an empirical analysis for a five-dimensional financial portfolio. Concluding remarks are given in 
Section 5. 



2 Smooth nonparametric vine copulas 

The purpose of this section is to shortly introduce the fundamentals of vine and Bernstein 
copulas. 



2.1 Vine copulas 

Vine copulas are obtained by hierarchical cascades of conditional bivariate copulas and are 
characterized by an increased flexibility for modeling inhomogeneous dependence structures in 
high dimensions. Here, we only present the basic definition as well as some basic properties of 
two popular classes of vines, i.e., C- and D-vinesE which will be used later on in our application to 
financial market data. Readers interested in a mo re rigorous exarn inati on of vine copulas and their 
prope rties are referred to the excellent studies by Joe ( 1996 . 1997 ) and Bedford and Cooke (jlOOl . 
2OO2I) . 

Starting point is the well-known observation that a joint probability density function of dimen- 
sion d can be decomposed into 



/(Xi, ...,Xd) = f{Xi) ■ f{x2\Xi) ■ /(XslXi, X2) ■ fiXd\Xi, . . . , Xd-l). 



(1) 



The classes of C- and D-vines are subsets of the so-called regular vines (or R-vines in short, see lBrechmann et al. 



2012h . We do not consider other types of R-vines in this paper but note that our proposed use of smooth Bernstein and 
B-spUne copulas as pair-copulas can also be extended to other subsets of R-vines. 
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Each factor in this product can then be decomposed further using a conditional copula, e.g. 

f{x2\x,) = Ci2(Fi(xi),F2(x2)) ■ /2(X2) (2) 

with Fj(-) being the cumulative distribution function (cdf) of Xi (i = 1, . . . , d) and Ci2(-) being the 
(in this case unconditional) copula density of (xi, X2). 

Going further down the initial decomposition, the conditional density /(xslzi, X2) could be 
factorised via 



/(X3|X1,X2) = C23|l(i^2|l(a;2|a;i),F3|i(x3|Xi)) ■ Ci3(Fi(xi), F3(X3)) ■ /3(X3) (3) 

with C23|i(-) being the conditional copula of (x2, a;3) given xi. 

Finally, substituting the elements of the initial decomposition in with the conditional copu- 
las, we get for dimension d = 3 

f{xi,X2,X3) = C23|i(F2|i(a;2|a;i),i^3|i(a;3|a;i)) (4) 

■ Ci2(Fi(xi),F2(x2)) 

■ Ci3(Fi(xi),F3(x3)) 

• fliXi) ■ /2(X2) ■ /3(X3) 

with C12, Ci3 and C23|i as pair-copulas. Note that as there are several possible decompositions of 
the conditional distributions, the joint density of Xi, . . . , a;^ can also be represented by different 



Bedford and Cooke 


( 


2001 




2002) 



dimensional joint density as a nested set of trees where two nodes are joined by an edge in tree 
j + l,j = l,...,(i— 1, only if the corresponding edges in tree j share a common node. Conse- 
quently, there are — 1 trees, where tree j has d + 1 — j nodes and d — j edges with each edge 
corresponding to a pair-copula density, i.e., a density of a conditional bivariate parametric copula. 
In a Canonical or C-vine, each tree has a unique node (without loss of generality this is node 1) 
that is connected to all other nodes yielding the representation 

d d~l d—j 

A;=l j=l i=l 

(5) 

where the subscript j identifies the tree, while i runs over all edges in each tree. In contrast to this. 
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in a D-vine, no node in any tree Tj is connected to more than two edges yielding the decomposition 



d d—1 d—j 

))• (6) 



k=l 



Examples of possible decompositions of a five-dimensional random vector via a C- and D-vine 
copula are shown in Figure [H and [2l respectively. 

— insert Figures [Hand |2] here — 

Fitting a vine copula model to a given dataset requires three separate steps. First, one needs to 
select the tree structure of the vine model. For a C- or D-vine, this amounts to the selection of a 
permutation of the indices 1, . . . , c/ of the random variables. As su ch, for a (j-dimen sional random 

Once a per- 



Aas et al 



200' 



vector their exist d\/2 different C- and D-vines, respectively (see 
mutation has been chosen, the structure of the vine is fully specified. Second, the statistician has to 
select d{d—l)/2 bivariate pair-copulas from candidate copula families. In the last step, the param 



DiBmann et al. 



eters of the pa ir-copulas have to be estimated. To select the optimal tree structure, 
dforthcomingl) propose a heuristic procedure in which the tree structure is chosen via a maximum 
spanning tree algorithm that maximizes the sum of the absolute empirical Kendall's r of all possi- 
ble variable pairs on a given level of the tree. Fo r the selection of the parametric pair-copulas. 



Brechmann and Czadd (120121) and 



Difimann et al, 



dforthcomingl) propose a sequential heuristic 



which selects the best fitting parametric copula family for each pair-copula based on the candi 
date copulas' AIC values. Although this sequential selection of param etric pair-copulas using 



Brechmann and Czado 



(120121) show that this heuristic yields considerably better results than a selection algorithm based 
on copula goodness-of-fit testsl 



AIC d oes not necessarily yield a globally optimal AIC value for the vine, 

f 

In the following, we concentrate on the problem of selecting the bivariate pair-copulas by 
substituting them with smooth nonparametric estimates of the underlying (pairwise) dependence 
structures. As a benchmark to ou r prop osed nonparametric method, we e mploy the sequential 



heuristic by Brechmann and Czadd (|2012h and 



Difimann et al. 



(|forthcoming| ). However, we expect 



our nonparametric approach to improve on the heuristic benchmark for several reasons. First, the 



nonparametric modeling of the pair-copulas eliminates comp 
incorrect parametric family for a given pair-copula. Second, 



etely the model risk o f choos ing an 



Grcinneberg and Hjorti (|2012r) prove 



^Choosing the best fittin g tree structure manually thus quickly becomes unfeasible in higher dimension (see 
DiBmann et al. .Iforthcoming). 



'^Results by iWeifi (2011) and iGrundke and Polld ( 120121) underline this finding. While IWeifil ( 120111) shows that 
goodness-of-fit tests give only little guidance for choosing co pulas when one is intereste d in forecasting quantiles in 
the extreme tails of a joint distribution, the empirical results of iGrundke and Polld (120 12h cast additional doubt on the 
ability of copula goodness-of-fit tests to identify stressed risk dependencies. 
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that the use of AIC as a model selection criterion is not correct in case rank-transformed pseudo- 



observations are used (as is common in almost af 
as already hinted at by 



Difimann et al. 



applications of copulas in finance). Third, 
(|forthcomingh . the incorrect specification of the parametric 



pair-copulas in the upper levels of a vine's tree structure can lead to a propagation and amplification 
of rounding errors causing the heuristic to become numerically unstable. 

In the next subsection, we define and discuss Bernstein copulas which we use as smooth non- 
parametric estimates of the pair-copulas in a vine model. 



2.2 Bernstein copulas 

As a nonparametric candidate for the pair-copulas in ^ and Q, we consider the recently 
proposed Bernstein copulas. In the following, we briefly state some basic mathematical facts on 
Bernstein polynomials and Bernstein copulas, respectively. The Bernstein polynomials of degree 
m are defined as 



5(m,A;,z) = (^^Jz*=(l -^)'"-^ (7) 

where k = 0, . . . ,m E N and < z < 1. 

As our focus lies on the nonparametric modeling of pair-copulas in vines, we restrict our anal- 
ysis in the following to bivariate Bernstein copulas. Let U = {Ui, U2) denote a discrete bivariate 
random vector with uniform margins over Tj := {0, 1, . . . rrii] with grid size mj G N and i = 1,2. 
In our analysis, we later choose nii := m2 := m = const. With 

2 

p{ku k^) := P{^{U, = h}), {ku k^) G [0, If (8) 

i=l 

we can then define the Bernstein copula density as 

mi — 1 m.2 — 1 2 

c(mi,M2) := ^ p{ki,k2)WmiB{mi-l,ki,Ui), (ui, M2) G [0, 1]^. (9) 

fci=0 fc2=0 j=l 



Pfeifer et al, 



(|2009|) show by integrating expression ^ that the cdf of the Bernstein copula is then 



given by 

C{xi,X2) c{ui,U2)duidu2 (10) 

Jo Jo 

mi m2 2 2 

ki=0k2=0 i=l i=l 

for (xi, X2) G [0, 1]^. Note that in order to smoothly approximate the distribution or density of a 
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copula in © and (fTOl) . very high degrees for the Bernstein polynomials have to be chosen. 

The Bernstein copula a s defined above ca n then be used to approximate the empirical copula 
process as defined, e.g., by lOeheuvelsl (|l979|) . To be precise, we approximate nonparametrically 
the joint distribution ofU = {Ui, U2) by using a bivariate sample {(Xj, Fj)}"^! of the underlying 
copula of s ize n . X(k) be the kih. order statistic of the sample. The empirical copula of 



Deheuvels 



(ll979Lll98lh is then defined as 



Cn{x,y) 




<y< 



(12) 



where = ^ ELi l{x,<X(„y,<y(,}, C„(04) = C„(i,0) = and C„(04) = 0(t,j = 

1,2, ...,72). 

To fit the Bernstein copula to the sample from the empirical copula, we first need to calculate 
the relative frequency of the observations in each target cell of a grid with given grid size m. The 
outcome of this is the contingency table [aki]k,i=i,...,m- Note, however, that the resulting marginals 
of the data [aki] do not need to be uniformly distributed and that the resulting approximation could 



Pfeifer et al. 



(120091) propose to transform 



therefore not be a copula. To circumvent this problem, 
the contingency table [aki] into a (possibly suboptimal) new contingency table [xki] with uniform 
marginals via a Lagrange optimization approach yielding 



a.j tti. 2 



— + ^ for i,j 



m 



m 



l,...,m. 



(13) 



where the index ■ denotes summation. Note that the quality of the Lagrange solution is reduced by 
an increasing number of the sample size n @We therefore chose to employ a different optimization 
strategy to correct for the non-uniform distribution of the marginals. 

Consequently, we calculate the approximation [xki] to the contingency table [aki] by solving 
the following optimization problem: 



subject to 



EE' 

k=l 1=1 



Xki — O-kl] 



mm 



(14) 



m m ^ 

^^Xfcj = ^^Xj; = — and > for j = 1, . . . , m. (15) 

k=l 1=1 ^ 



To solve for the \xk i], we make use of the quadratic optimization algorithm of 

Goldfarb and Idnanil (|l982l . Il983r) . In preliminary tests, the found solutions to this optimiza- 



In unreported results, the optimization strategy of lPfeifer et al.l (120091) proved to yield only suboptimal results. 
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tio n problem y i elded significantly lower quadratic errors than the procedure initially proposed 
by 



Pfeifer et al, 



(|2009l) thus confirming the need for a more refined optimization strategy. 
To use Bernstein copulas as pair-copulas both in our simulation study and the empirical appli- 
cation, we require efficient algorithms for simulating and evaluating the density a nd distribution of 

J2OO9I) 



Aas et al 



a given vine copula. To this end, we adapt the algorithms initially proposed by 
by substituting the parametric h-hunctions (i.e., the partial derivatives of the copula densities) in 
these algorithms by the partial derivatives of the fitted bivariate Bernstein copulas. 



3 Simulations 



In this section, we illustrate the superiority of t he smooth nonparametric vin e model over the 



sequential heuristic of 



Brechmann et al 



Difimann et al 



(|forthcomingl) for selecting the 



mm and 

pair-copulas in a vine parametrically. In particular, we are interested in the error of the approxima- 
tions to a pre-specified true copula using both our smooth nonparametric model as well as a vine 
model calibrated \ yith parametric pa ir-copulas. The setup of our simulation study follows the pro- 



cedure laid out in 



Shen et al. 



(120081) . but differs in that way that we also consider a (heuristically 
calibrated) parametric benchmark approximation to the true vine copula. 

As a measure for the approximation error, we compare the pre-specified true pair-copulas of 
the vine model with the parametric and nonparametric approximations and use the average squared 
error (ASE) of the cdfs of all bivariate pair-copulas each taken at mi x 7712 uniform grid points in 



, I.e., 



ASE :-- 



d{d — 1) nil ■ 



d{d-l)/2 mi m2 

E EE (6. 

i=i j=i k=i 



J 



k 



nil + 1 1712 + 1 



J 



k 



nil + 1 + 1 



(16) 

where Cj is the cdf of a pre-specified pair-copula from which we simulate a random sample of size 
n and Ci is a (parametric or nonparametric) approximation to the pair-copula Ci computed on the 
basis of the simulated random sample. 

In the simulation study, we consider two sample sizes n = 200 and n = 500 to assess the de- 
creasing effect of the sample size on the approximation error. Furthermore, we analyze the effect of 
the type (C- or D-vine) as well as the dimensionality of the vine model on the approximation errors. 
To be precise, we simulate random samples from vines of dimension d = 3, 5, 6, 7, 11, 13. As the 
dimension of the vine model increases, so does the number of variables one has to condition on in 
the pair-copulas of the vine's lower trees. The pair-copulas in the lower trees of the vine, however, 
are generally more complex to model so that the accurate approximation of the pair-copulas on all 
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levels of the vine constitutes a considerable challenge to our nonparametric approximation At 
the same time, the curse of dimensionality could additionally complicate the approximation of the 
pair-copulas thus making the comparison of our approximation for different dimensions a sensible 
exercise. Finally, we expect the propagation and amplification of rounding errors to increase in 
higher dimensions possibly leading to the numerical instability of the parametric heuristic. 

As candidate parametric copula families from which the pair-copulas of the true vine models 
are chosen, we use the Gaussian, Student's t, Clayton, Gumbel, Survival Clayton, Survival Gum- 
bel, the rotated Clayton copula (90 degrees) and the rotated Gumbel copula (90 degrees). The 
pair-copulas as well as their respective parameters in each simulation are chosen randomly. For 
each sample size, dimension and vine type, we simulate 1, 000 random samples and approximate 
the data with a vine copula using Bernstein copulas as pair-copulas. As we are only interested in 
measuring the accuracy of the approximation of the pair-copulas, we calibrate our nonparametric 
vine using the correct vine type as well as the correct tree structure. As a benchmark, we call 
brate a sec ond vine copula by using the sequential heuristic proposed by 



(2012) and 



Difimann et al. 



Brechmann and Czado 



(IforthcomingJ) . Furthermore, we also compute the fraction of time the 
sequential procedure breaks down due to either the numerical instability of the evaluation of the 
likelihood function and the parameter estimation or due to the ASE tending to infinity. The results 
of the simulations are presented in Tabled] 

— insert Tabled] here — 

The results shown in Table [I] present several interesting insights into the finite sample prop- 
erties of both the heuristically calibrated parametric and our proposed nonparametric vine copula 
models. First, we can see from Tabledlthat for lower dimensions (e.g., d = 3 and d = 5) the ASE 
of our nonparametric is considerably larger than for the parametric model calibrated by sequen- 
tially selecting the pair-copulas based on AIC values. With increasing dimension of the random 
vector, however, we can observe that the approximation error of the parametric model increases 
disproportionately compared to our proposed nonparametric model. Furthermore, the nonpara- 
metric model appears to be able to match the approximation error of the parametric approach for 
dimensions d = 13 and higher. Most importantly, the parametric modeling approach becomes 
highly numerically unstable in higher dimensions. At the same time, our proposed nonparametric 
vine with Bernstein pair-copulas is extremely reliable yielding acceptable approximations to the 
true underlying dependence structure even for high-dimensional random vectors. The parametric 



"'This is one reason whv lAas et al. ( 2009b . Brechmann et al.l ( 2012 ) and DiBmann et al. ( forthcomind) propose to 



capture as much dependence of the joint distribution that is to be modeled in the first trees of a vine model. If these 
pair-copulas are modeled accurately, the remaining pair-copulas in the lower trees can then be truncated or simplified. 
Furthermore, the truncation and simplification of a vine on the lower levels of the vine's tree limits the potential 
propagation of rounding errors. 
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approach, on the other hand, breaks down in approximately 50% of all simulations for dimension 
d = 13 and higher. In many of these cases, the bad approximation (or numerical instability) of 
the parametric approach was caused by the wrong selection of several parametric families for the 
pair-copulas in the vine model. Concerning the type of the vine copula model, we find no signifi- 
cant differences between the average approximation errors of the C- or D-vines. As expected, we 
also find the average approximation error of both the parametric and nonparametric model to be 
decreasing in the sample size used for estimating both models. 

To further illustrate the finding that the nonparametric model improves on the accuracy of a 
parametric vine especially in higher dimensions, we plot simulated samples from several paramet- 
rically and nonparametrically fitted copulas in Figure [3] where we assume that the true underlying 
dependence structure is given by a Clayton copula with parameter 9 = 5. From this copula, we 
simulate a random sample of size n = 500 and fit both a nonparametric Bernstein copula as well 
as a parametric Clayton and Gumbel copula via Maximum-Likelihood to the data. From all three 
fitted copulas, we again simulate a random sample and compare the plots of the simulated obser- 
vations with the original sample. 

— insert Figure [3] here — 

The plots in Panels (a) and (b) in Figure [3] show how the Bernstein copula approximately cap- 
tures the lower tail dependence of the original sample shown in Panel (a). The plot of the Bernstein 
copula in Panel (b) also shows, however, that the nonparametric approximation of the data sam- 
ple coincides with a loss in information on the tail behaviour of the true underlying dependence 
structure. At the same time. Panel (c) underlines the notion that the nonparametric model is not 
superior to a parametric model in which the parametric copula family has been chosen correctly. 
If, however, the parametric copula familiy is chosen incorrectly like it is shown in Panel (d), the 
wrong selection of the parametric copula family can cause considerable approximation errors and 
a severely inaccurate modeling of the underlying tail dependence. Given no prior information on 
the parametric copula familiy, the nonparametric Bernstein copula model clearly improves on the 
fit of an inaccurately fitted parametric model. The nonparametric modeling of the pair-copulas 
thus seems to be a sensible approach especially when the number of pair-copulas that need to be 
selected from candidate parametric copula families increases (i.e., with increasing dimension). 

4 Empirical study 
4.1 Methodology 

The purpose of our empirical study is to investigate the superiority of the proposed smooth 
nonparametric vine copula models over the competing calibration strategy based on a sequential 
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selection of parametric pair-copulas via AIC with regard to the accurate forecasting of a portfo- 
lio's VaR. The simulation study presented in the previous section has highlighted the finding that 
our vine model with smooth nonparametric pair-copulas is especially well-suited for dependence 
modeling in higher dimensions as the selection of parametric pair-copulas becomes numerically 
unstable due to error propagation and amplification. For low-dimensional problems, the heuristic 
selection of parametric copulas, however, seems to outperform our nonparametric approach with 
respect to the ASE of the vine copula's approximation. To show that our nonparametric model 
matches the results of the parametric heuristic even for low-dimensional problems, we concentrate 
in our empirical analysis on the VaR forecasts of a five-dimensional portfolio. 

Financial data are usually characterized by the presence of both conditional heteroscedas- 
ticity and asymmetric dependence in the log returns on financial retu rns. Therefore, we fol- 
low the vast mai ority of studies on copula models for VaR-estimation ( Jondeau and Rockingei . 



20061 : iFantazzinil |2009|; lAusin and Lopesl l2010l : iHafner and Reznikoval 



20101) and employ stan- 



dard GARCH(l,l)-models with Student's t-distributed innovations to model the marginal be- 
haviour of our data. Although different specifications of the GARCH model are also possible, 
results found by lHansen and Lundd (|2005|) suggest that the choice of the order of a GARCH model 
is only of little importance for the model's forecasting accuracy. 

Throughout the empirical study, we consider continuous log returns on financial assets with 
prices Pt {t = 0, 1, . . . ,T). The assets' log returns Rt are defined by Rt := \og{Pt/ Pt-i) for 
t > 1. Our focus lies on modeling the joint distribution of the d assets, i.e., the joint distribution of 
the returns Rn, . . . , Rtd- 

The marginal behaviour of the assets is modeled by the use of GARCH(l,l)-models with t- 
distributed innovations. The marginal model is then given by 



R, 



fXj + atjZtj, 

aoj + aijRf_^j + j = l,...,d; t = l,...,T, 



(17) 
(18) 



with independent and identically t-distributed innovations Ztj. The dependence structure between 
the d assets is introduced into the model by assuming the vector Zt = {Za, . . . , Ztd) (t = 1, . . . , T) 
of the innovations to be jointly distibuted under a d-dimensional copula C with 



(19) 



where i^i, . . . , i/^ are the parameter vectors of the innovations, C is a copula and a; is a vector of 
copula parameters (in case of the parametric model, otherwise is simply empty). 

The parameters of the univariate GARCH-models are estimated via Quasi-Maximum Likeli- 
hood Estimation. For the estimation of both the nonparametric vine model as well as the para- 
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metric model calibrated by using the pair-copulas' AIC values, we make use of rank-transformed 
pseudo-observations rather than the original sample as input data0 As the main results for copu- 
las only hold for i.i.d. samples, we use the parameter estimates for the univariate GARCH models 
and transform the original observations into standardize d residuals to yield (approx imately) i.i.d. 
observations before computing the pseudo-observations (|Dias and Embrechts . 



2009 ). 



In our empirical application, we consider an equally-weighted five-dimensional portfolio with 
returns Rp t = X]j=i ^tj- The results from our simulation study underline the finding that our 
proposed vine copula model with Bernstein pair-copulas, on average, yields a better approximation 
to the empirical copula than the heuristically calibrated parametric model especially in higher 
dimensions. However, the parametric copula vine model could still outperform our proposed model 
w.r.t. the forecasting accuracy in low dimensions. In our empirical application, we therefore restrict 
our analysis to a portfolio consisting of five assets to additionally illustrate the nonparametric 
Bernstein vine copula model's superiority for low-dimensional problems. 

To forecast the portfoli o returns, we employ the algorithm presented in the study by 
Nikoloulopoulos et al.l (120111) initially pro posed for in-sample forecasting which was extended to 



WeiJ3l ( 2012 ) 



out-of-sample forecasting by 

The aim of the algorithm is the computation of a one-day-ahead forecast for the portfolio re- 
turn i?p f via Monte Carlo simulation. In a first step, K = 10, 000 observations • • • , ^t+i d 
(k = 1, . . . , K) from the fitted (parametric or nonparametric) vine copula are simulated. Us- 
ing the quantile function of the fitted marginal Student's t distributions, the simulated vine cop- 
ula observations are then transformed into observations z^^-^ ■ from the joint distribution of the 
innovations. In the next step, the simulated innovations are transformed into simulated returns 



-Rr+1 J = + ^T+ijZT+ij where (3"t+i j and fij are the forecasted conditional volatility and mean 
values from the previously fitted marginal GARCH models. The MC-simulated forecasts of the 
portfolio return is then simply given by -R^+i ^ = Yl'^j=i ^t+ij- Sorting the simulated portfo- 
lio returns for a given day in the forecasting period and taking the empirical one-day a percentile 
then yields the forecasted a%-VaR. 

To backtest th e results of our for ecasting, we employ the test of conditional cover- 
age proposed by IChristoffersenI (Il998h as well as two duration-based tests discussed in 
Christoffersen and Pelletiei ( 2004a ). ^ 



"For a comp arative study on the finite sample properties of different ML-based estimators for copulas, see 
Kim et al.l ( 120071) . The authors show that absent any information on the true distribution of the marginals, statisti- 



cal infer ences should be based o n rank-transformed pseudo-observations. 

'^See lBerkowitz et al.l (1201 II) for an excellent review of different me thods for bac ktestin g Valu e-at-Risk forecasts. 
A comparison of different backtests can be found in the recent study by lEscanciano and Peii (120121) . 
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All three backtests are based on the hit sequence of VaR-exceedances which is defined by 



1, if <VaR„ (i?p,t) I 



0, otherwise 



with t being the ti me subscript and bein g the set of available inforrnation. T he test of condi 
tional coverage by 



ChristoffersenI (11998 ) and 



Christoffersen and Pelletierl (|2004al) jointly tests for 



the correct number of VaR-exceedances (unconditional coverage) and the serial independence of 
the violations over the complete out-of-sample (independence)!^ Under the null hypothesis of 
a correct num ber of VaR-exceedances that are ind ependent over time, the hit sequence is simply 
distributed as ( Christoffersen and Pelletierll2004a ) 



''t,a 



i.i.d. Bernoulli{a). 



Then, let P be the length of the out-of-sample. Pi be the number of VaR-exceedances and Pq be the 
number of days on which the daily VaR forecast was not exceeded, respectively (and consequently 
P = Pi + Pq). The likelihood function for the i.i.d. Bernoulli hit sequence with unknown 
parameter tti is 

L{ht,o.,TTi) = n^'{l-7iif-''^ (20) 

and the Maximum-Likelihood estimate of tti is simply given by tti = Pi/ P. The test of uncondi- 
tional coverage is then given by a likelihood ratio test based on 



LRuc = -2 (In L(/ij,„, tti) - In L{ht,a, a)) . 



(21) 



To test the hypothesis of independently distributed hits, the hit sequence is assumed to follow a 
first order Markov sequence with switching probability matrix 



n 



1 - TToi TToi 
1 - TTii TTii 



(22) 



with TTjj being the probability of an i on day t — 1 being followed by a j on the next day t and 
i,j E {1; 0}. Using the likelihood function 



■JT, 



01 



TTll) 



Pi-Pii^Pi 
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'^The test of unconditional covera ge has been implicitly incorporated in the Basel A ccord for determining cap- 
ital requirements for market risks, see iBasel Committee on Ban king Supervision' 

(■1994. 

Consequently, it has since 
become an industry standard in market risk management, see, e.g.. iEscanciano and Peii (I2012h . 
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where Pij is the number of observations in hf^a where a j follows an i and i,j G {1; 0} and the 
ML-estimates ttqi = P01/-P0 and itn = Pn/Pi, the likelihood ratio test of the independence of 
hits is given by 

LRind = 2 (lnL(/ij,o,7roi,^ii) - \nL{ht,aAi)) ■ (23) 

Both tests are then combined via LRcc = LRuc + LRind to yield the test of conditional cov- 
erage. We note here that we do not rely on the asymptotic C hi-squared distribution of the test 



HsuetaL 



(120 lib . Although easy to implement, p- 



statistic which is used, e.g., in the study by 
values derived under the assumption of the test statistic following a Chi-squared distribution are 
usuall y incorrect due to the gener a lly low sample sizes when using hit sequences. Instead, we 
follow IChristoffersen and Pelletieii (|2004ah in generating approximate p-values via Monte Carlo- 



simulation. 



As an alternative to the test of conditional coverage, IChristoffersen and Pelletien (|2004a|) pro- 
pose backtests based on the durations between VaR-exceedances. Then, let 



U-t 



i-l 



(24) 



Christoffersen and Pelletierl (l2004al) propose to use 



^bD^-'^ exp (-(aD)^) which nests 



1. 



be the duration of time (in trading days) between two subsequent VaR-exceedances where ti is the 
time of the ith VaR-exceedance. Under the null hypothesis of a correctly specified VaR model, 
we would expect the process of no-hit durations to have no memory and mean 1/a. Conse- 
quently, the process D of durations should fo llow an exponential distribu t ion wit h fexpiD', ol 
a exp (— a;D)o As an alternative hypothesis, 
the WeibuU distribution for the process D with fw{D; a, b) 
the exponential distribution from the null hypothesis for b - 

Although this test potentially captures higher order dependence in the hit sequence /i^ „, the 
informati on from the temporal ordering of the no-hit durations is not exploited in the backtest. As 
a remedy, IChristoffersen and Pelletien (|2004ar) propose a conditional d uration-based test based o n 
the Exponential Autoregressive Conditional Duration (EACD) model of lEngle and Russelll (| 19981) . 
In the standard EACD (1,0) model, the conditional expected duration Ej_i(Dj) is assumed to 
follow the process 

Ei„i(A) = W + (25) 
Again assuming an underlying exponential density with mean equal to one in the null hypothesis. 
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See Christoffersen and Pelletieii (l2004al) for details of the backtest and the motivation for using a continuous 
distribution for the discrete process D. 
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the conditional distribution of the duration is is given by 

Ieacd (Dilipi) = ^ exp (^--^^ ■ (26) 
The null hypothesis of independent no-hit durations is then given hy Hq : (3 = 0. 

4.2 Data 

We consider a five-dimensional equal- weighted portfolio consisting of the returns on the EURO 
STOXX 50 Price Index, 30-year US Treasury Bonds, France Benchmark 10-year Government 
Bond Index, Gold Bullion LBM and one-month forward Crude Oil Brent. We obtain the data 
from the Thomson Reu ters Datastream database. We follow the screening procedure proposed by 



Ince and Porteii (120061) to control for known sources of data errors in Datastream. To be precise, 
we check whether our data include prices below $ 1 (which could lead to erroneous log returns due 
to Datastream'?, practice of rounding prices) as well as log returns above 300% that are reversed 
within one month. Our five univariate time series do not suffer from any of these data errors. 

Our sample covers a period of 800 trading days ranging from June 6, 2009 to July 19, 2012 and 
thus includes the aftermath of the default of Lehman Bros, as well as the onset of the Sovereign 
Debt Crisis. We use rolling windows with a length of 500 trading days for forecasting the one- 
day-ahead VaR on the following trading day. Our full out-of-sample spans a period of 300 trading 
days. Time series plots of the five portfolio constituents as well as the returns on the equal-weighted 
portfolio are shown in Figure IH Panels (a) through (e) show the time series plots of the univariate 
returns, while Panel (f) shows the time series plot of the portfolio. The initial in-sample is shaded 
in grey to highlight the out-of-sample consisting of 300 trading days. 

— insert Figure |4] here — 

The plots in Figure|4]show several distinct features that can complicate VaR- forecasting. First, 
all plots exhibit the common stylized fact of volatility clusters, e.g., in Panels (a) and (c). Second, 
overall volatility of the univariate returns differs significantly across our five portfolio constituents. 
For example, while the returns on the 30-year US treasury bonds are quite volatile in the in-sample 
and seem to calm in the out-sample, the opposite is true for the France Benchmark 10-year Govern- 
ment Bond Index which exhibits low volatility in the in- sample and a pronounced cluster of high 
volatility and extreme spikes around November 2011. This last result is, however, not surprising 
considering the fact that the Sovereign Debt Crisis experienced a climax at that time with the res- 
ignation of the Greek and Italian Prime Ministers, early elections in Spain and the expansion of the 
European Financial Stabilisation Mechanism (EFSM). Similarly, the price of Gold bullion became 
more volatile in the out-of-sample as well. The plot in Panel (f) shows that the combination of 
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the five individual assets produces a portfolio which exhibits several phases of both high and low 
volatility as well as sudden extreme spikes in the portfolio's log returns. 

4.3 Results 

Following the methodology presented in Section |4~T1 we compute the one- step-ahead forecasts 
of the portfolio- VaR for each day in the out-of-sample using rolling windows of 500 trading days. 
To analyze the differential effect of different confidence levels for the VaR on our models' fore- 
casting accuracy, we forecast the 2%-, 5%- and 10%-VaR for a long and the 97.5%-VaR for a short 
position in the portfolio. Thus, we would expect 6, 15, 30 and 8 exceedances below the forecasted 
VaRs, respectively^ The VaR-forecasts as well as the realized portfolio returns for all four con- 
fidence levels are shown in Figures |4] and [5l In both figures. Panels (a) and (b) show the realized 
portfolio returns and the VaR forecasts computed by the use of the nonparametric vine copula 
and the parametric vine copula model calibrated via the heuristic based on Akaike's Information 
Criterion, respectively. 

— insert Figures |4] and [5] here — 

The plots in Figures |4] and [5] show that both the nonparametric and the parametric model yield 
rather accurate forecasts of the portfolio's risk. While all VaR-forecasts are sufficiently close to 
the realized portfolio returns, exceedances of the VaR-forecasts occur only in case of large losses 
on the portfolio investment. Also, we can see that both the parametric and nonparametric model 
specifications yield quite similar VaR-forecasts. Thus, the proposed nonparametric vine copula 
model seems to perform exceptionally well even for low-dimensional portfolios. Furthermore, 
the finding of comparable VaR forecasts of both the nonparametric and parametric model remains 
valid for all four VaR confidence levels we consider. In addition to this, the results presented in 
Figures |4] and [5] also show that both models adequately adapt the VaR-forecasts to changes in the 
portfolio returns' volatility^ To further assist in the interpretation of the results. Figures [6] and |7] 
highlight the positive and negative VaR-exceedances for all models and the three confidence levels 
for a long position in the portfolio. 

— insert Figures [6] and |7] here — 

Both figures underline our first impression from Figures |4] and [5] that both models forecast 
the VaR of the portfolio quite accurately. We can see from Figure \5\ that not only do both mod- 
els yield (approximately) correct numbers of VaR-exceedances for all three confidence levels, the 

'^For the short position, VaR-exceedances are defined as returns above the daily VaR forecast. 

'^The flexible adjustment of both models to changes in return volatility also underlines the fact that a static de- 
pendence model in conjunction with dynamic marginal models suffices to model and forecast the dynamics of a 
multivariate return distribution. 
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exceedances also seem to occur randomly in time. Most importantly, however, our proposed non- 
parametric vine copula model with GARCH margins easily matches the heuristically calibrated 
parametric vine w.r.t. the accuracy of VaR-forecasting even for a relatively low-dimensional port- 
folio. To further substantiate this finding, we perform three formal backtests on the results of both 
the parametric and nonparametric vine models. The results of the three backtests are presented in 
Tabled 

— insert Table [21 here — 

The backtesting results stress our finding that both models yield comparable results. For ex- 
ample, all but one VaR models cannot be rejected at the 99% confidence level based on the test of 
conditional coverage. Although the results of the unconditional duration-based backtest imply a 
significantly worse forecasting accuracy of both models, the p-values for both the nonparametric 
and the parametric model are comparable for different confidence levels of the VaR. This indicates 
that neither model outperfo rms the other one based on our seco nd backtest. If we use the condi- 
tional duration-based test of IChristoffersen and Pelletien (|2004bf) instead, none of the VaR-models 



is rejected. Turning to the number of VaR-exceedances, the results of our nonparametric vine cop- 
ula model are slightly better for the 5%-VaR than those of the parametric benchmark while the 
opposite is true for the (for most practical uses too optimistic and thus unsuitable) 10%- VaR. 

Our backtesting results indicate that both models yield acceptable VaR- forecasts for a relatively 
low-dimensional portfolio. One could conclude from this finding that in general using our non- 
parametric vine copula model does not yield significantly better VaR- forecasts. However, one has 
to keep in mind that our empirical analysis was deliberately aimed at testing the hypothesis that the 
nonparametric model yields accurate VaR-forecasts even in lower dimensions. In unreported tests 
of high-dimensional portfolios, the parametric benchmark suffered from the same numerical insta- 
bility that was also observed in our simulation study. At the same time, our nonparametric model 
produced accurate VaR-forecasts in a numerically stable fashion even for high-dimensional port- 
folios when the parametric benchmark had either broken down or produced woefully inaccurate 
VaR-forecasts. 



5 Summary 

In this paper, we propose to model the pair-copulas in a vine copula model nonparametrically 
by the use of Bernstein copulas. Our proposed model has the advantage of a significantly reduced 
model risk as it avoids the error-prone selection of pair-copulas from candidate parametric copula 
families. In contrast to previous studies on the use of discrete empirical copulas as pair-copulas, our 
proposed use of Bernstein copulas has the additional advantage that the building blocks in a vine 
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model are approximated by smooth functions from which one can easily simulate random samples. 
We test the approximation error of the smooth nonparametric Bernstein vine copula model against a 
parametric benchmark calibrated by the use of a sequential heuristic based on AIC. The superiority 
of our proposed model is exemplified in an empirical risk management application. 

The results we find in our simulation study show that for low-dimensional problems, the para- 
metric modeling approach outperforms our proposed nonparametric approach only marginally. 
However, the differences in the approximation error quickly vanish for higher dimensions with 
both models yielding comparable average approximations errors for dimensions d = 13 and higher. 
At the same time, our proposed nonparametric vine copula model does not suffer from numerical 
instability and error propagation which plagues the parametric benchmark due to an increasing 
number of wrongly selected parametric pair-copulas. 

In the empirical risk management application, we test whether the differences in the average 
approximation error of the parametric and nonparametric vine copula models cause significant dif- 
ferences in both models' accuracy of forecasting the VaR of a low-dimensional asset portfolio. The 
results of our analysis show that even in lower dimensions (d = 5), our nonparametric vine cop- 
ula model yields VaR- forecasts that cannot be rejected by several different formal backtests. The 
proposed nonparametric vine copula model thus seems to match the (good) results of a parametric 
vine copula model in lower dimensions and significantly outperforms this benchmark in higher 
dimensions. 

A natural extension of our model would be to consider more sophisticated smooth approxi- 
mations of the empirical copula. Cubic B-splines and non-uniform rational B-splines (NURBS) 
appear as natural candidates for this job. While Bernstein copulas have been shown to be good 
smooth nonparametric replacements for parametric pair-copulas, spline copulas should yield even 
better approximations while at the same time yielding numerically stable vine model calibrations 
as well. We intend to analyze the suitability of spline copulas in vine models in future research. 
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Figure 1: Five-dimensional C-vine copula. The figure shows an example of a five-dimensional C-vine 
copula with five random variables, four trees and ten edges. The nodes in the first tree correspond to 
the five random variables that are being modeled and each edge corresponds to a bivariate conditional or 
unconditional pair-copula. 



24 




Figure 2: Five-dimensional D-vine copula. The figure shows an example of a five-dimensional D-vine 
copula with five random variables, four trees and ten edges. The nodes in the first tree correspond to 
the five random variables that are being modeled and each edge corresponds to a bivariate conditional or 
unconditional pair-copula. 
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Figure 3: Simulated observations from parametrically and nonparametrically fitted copulas. The fig- 
ure shows plots of n = 1000 observations simulated from different parametric and nonparametric copulas. 
Panel (a) shows the plot of the original observations simulated from a Clayton copula with parameter 9 = 5 
which are used to calibrate a nonparametric Bernstein copula and two parametric Clayton and Gumbel copu- 
las. Panel (b) shows a sample of simulated observations from a Bernstein copula which was calibrated based 
on the sample in Panel (a). Panels (c) and (d) show similar plots of simulated samples from a parametric 
Clayton and Gumbel copula which were fitted via Maximum-Likelihood using the original sample shown 
in Panel (a). 
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Figure 4: Time series plots of log returns in the full sample. The figure shows plots of the log returns on the EURO STOXX 50 Price Index in 
Panel (a), US Treasury Bonds (30-year) in Panel (b), France Benchmark 10- Year Government Bond Index (Clean Price) in Panel (c), Gold Bullion 
LBM ($/Troy Ounce) in Panel (d), Crude Oil-Brent one-month forward ($/BBL) in Panel (e) and the returns on an equal-weighted portfolio consisting 
of the five individual assets in Panel (f). The sample covers the period from June 15, 2009 to July 19, 2012 (800 ti'ading days). The plots show the 
log returns during our complete sample and are divided into the initial in-sample of 500 trading days (shaded in grey) and the out-of-sample of 300 
trading days. 
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Actual portfolio returns and VaR forecasts: 
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Figure 4: Comparison of 10%— and 97.5%— VaR forecasts and realized portfolio returns. The figure shows plots of the log returns on the 
five-dimensional portfolio we consider in our empirical study and the VaR forecasts. Panel (a) presents the realized portfolio returns and the VaR 
forecasts computed by the use of the nonpaiametric vine copula as the dependence model. Panel (b) shows a corresponding comparison of the 
portfolio returns and the VaR forecasts estimated via a vine copula with the parametric pair-copulas chosen according to the sequential heuristic 
taken from the R-package CDVine based on Akaike's Information Criterion. For both models, the (1 — a)-VaR is computed for confidence levels 
a € {10%; 97.5%}. Both plots show results for the out-of-sample of 300 trading days. The portfolio consists of the returns on the EURO STOXX 50 
Price Index, US Treasury Bonds (30-year), France Benchmark 10- Year Government Bond Index, Gold Bullion LBM and Crude Oil-Brent one-month 
forward. 
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Figure 5: Comparison of 2%— and 5%— VaR forecasts and realized portfolio returns. The figure shows plots of the log returns on the five- 
dimensional portfolio we consider in our empirical study and the VaR forecasts. Panel (a) presents the realized portfolio returns and the VaR forecasts 
computed by the use of the nonparametric vine copula as the dependence model. Panel (b) shows a corresponding comparison of the portfolio returns 
and the VaR forecasts estimated via a vine copula with the parametric pair-copulas chosen according to the sequential heuristic taken from the R- 
package CDVine based on Akaike's Information Criterion. For both models, the (1 — a)-VaR is computed for confidence levels a G {2%; 5%}. Both 
plots show results for the out-of-sample of 300 trading days. The portfolio consists of the returns on the EURO STOXX 50 Price Index, US Treasury 
Bonds (30-year), France Benchmark 10- Year Government Bond Index, Gold Bullion LBM and Crude Oil-Brent one-month forward. 
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Absolute VaR (2%) Exceedances: 
(a) Bernstein-Vine 
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Figure 7: Negative VaR-exceedances for the Bernstein Vine and the parametric benchmark. The figure shows plots of the negative VaR- 
exceedances (i.e., losses below the daily VaR-forecasts) computed from the nonparametric Bernstein vine copula model (Panels (a), (c) and (e)) 
and the parametric benchmark vine model with the parametric pair-copulas chosen according to the sequential heuristic taken from the R-package 
CDVine based on Akaike's Information Criterion (Panels (b), (d) and (f)). For both models, the (1 — a)-VaR is computed for confidence levels 
a S {2%; 5%; 10%}. Both plots show results for the out-of-sample of 300 trading days. The portfolio consists of the returns on the EURO STOXX 
50 Price Index, US Treasury Bonds (30-year), France Benchmark 10- Year Government Bond Index, Gold Bullion LBM and Crude Oil-Brent one- 
month forward. For ease of presentation the losses exceeding the VaR-forecasts ai^e shown as positive real numbers. 



Table 1: Comparison of the Average Squared Error. The table presents a comparison of the Average 
Squared Errors (ASE) of the parametric and nonparametric approximation to randomly specified vine copu- 
las as well as the fraction of time (in %) the sequential heuristic based on AIC breaks down due to either the 
numerical instability of the evaluation of the likelihood function and the parameter estimation or due to the 
ASE tending to infinity. The ASE is given in multiples of 10~^. All results are given in averages of 1, 000 
simulations. 
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C-Vine {n = 
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0.040336 


6.90 


2.741199 
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0.030566 
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0.00 


Dimension d 


= 9 










C-Vine {n = 


200) 


5.394219 


37.40 


6.816040 
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Table 2: Backtesting results. The table presents the resuks of three different backtests performed on the out-of-sample forecasts for the portfoho 
VaR estimated from the vine copula models calibrated pai^ametrically via the sequential selection of pair-copulas via AIC and the nonpai^ametric 
modeling of the pair-copulas using Bernstein c opulas, respectively. The three back tests are the test of conditional coverage, the unconditional and 
the conditional duration-based tests proposed in lChristoffersen and Pelletierl (l2004al) . The table reports the expected and the realized number of VaR- 
exceedances as well as the p-values for the three backtests. For both models, the backtesting results are reported for the (1 — a)-VaR for confidence 
levels a € {2%; 5%; 10%; 97.5%}. For the 97.5%-VaR, exceedances are given under the assumption of a short position in the portfolio. 



Sequential AIC Bernstein Pair-Copulas 
Exceedances P-value Exceedances P-value Exceedances 
(expected) (realized) (realized) 

VaR a = 2% 
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Unconditional Duration 
Conditional Duration 
VaR a = 5% 
Conditional Coverage 
Unconditional Duration 
Conditional Duration 
VaR a = 10% 
Conditional Coverage 
Unconditional Duration 
Conditional Duration 
VaR a = 97.5% 

Conditional Coverage 8 0.116 3 0.117 3 

Unconditional Duration 8 0.053 3 0.052 3 

Conditional Duration 8 1.000 3 1.000 3 



6 0.994 6 0.850 5 

6 0.017 6 0.011 5 

6 0.999 6 0.999 5 



15 0.437 16 0.573 15 

15 0.424 16 0.186 15 

15 0.166 16 0.132 15 



30 0.152 34 0.042 37 

30 0.028 34 0.039 37 

30 0.991 34 0.990 37 



